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Abstract: The recovery of a spherically-symmetric wave speed v is considered in a 
bounded spherical region of radius b from the set of the corresponding transmission eigen- 
values for which the corresponding eigenfunctions are also spherically symmetric. If the 
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integral of 1/v on the interval [0, b] is less than 6, assuming that there exists at least one 
v corresponding to the data, it is shown that v is uniquely determined by the data con- 
sisting of such transmission eigenvalues and their "multiplicities," where the "multiplicity" 
is defined as the multiplicity of the transmission eigenvalue as a zero of a key quantity. 
When that integral is equal to 6, the unique recovery is obtained when the data contains 
one additional piece of information. Some similar results are presented for the unique 
determination of the potential from the transmission eigenvalues with "multiplicities" for 
a related Schrodinger equation. 
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1. INTRODUCTION 

The interior transmission problem is a nonselfadjoint boundary- value problem for a 

pair of fields \1/ and in a bounded and simply connected domain O of R n with the 

sufficiently smooth boundary dfl. It can be formulated as 

' A^ + Ap(x)^ = 0, xefi, 

on on 

where A denotes the Laplacian, A is the spectral parameter, n represents the outward unit 
normal to the boundary dfl, and the positive quantity p(x) corresponds to the square of the 
refractive index of the medium at location x in the electromagnetic case or the reciprocal 



of the square of the sound speed v(x) in the acoustic case, i.e. u(x) := 1/ sj /o(x). In the 



acoustic case, a/p(x) is usually called the slowness. Without loss of generality we can 
assume that in the region exterior to fi, the speed of the electromagnetic wave is 1 or the 
sound speed is 1 in the acoustic case. 

This interior transmission problem arises in the inverse scattering theory in inhomoge- 
neous media, where the goal is to determine the function p in O from an appropriate set of 
A- values related to (1.1). The values of A for which (1.1) has a pair of nontrivial solutions 
\1/ and \l/o are called transmission eigenvalues. It is already known that those transmission 
eigenvalues can be determined from some far-field measurements (see e.g. [5,7,13-15,31]). 

Since there is not a standard theory to analyze nonselfadjoint eigenvalue problems, the 
existence of transmission eigenvalues for (1.1) was an open problem until recently. Using 
some techniques related to the Fredholm theory of integral equations, it has been shown 
[14] that the transmission eigenvalues for (1.1) form a discrete set with infinity as the only 
possible accumulation point. In general we expect transmission eigenvalues to be complex 
numbers although some of them may be real and some, in fact, may be positive. 

Under the assumption that p(x) > 1 for all x e O (or the assumption that p(x) < 

2 



1 for all x G O), it has been shown [11] that the corresponding positive transmission 
eigenvalues for (1.1) form a countably infinite set. Similar results for the existence of a 
countably infinite set of transmission eigenvalues have been obtained [9,10,12,18,28] for 
related problems involving Helmholtz and Maxwell's equations, where the bounded region 
O is allowed to contain cavities, where p(x) = 1 in each cavity. 

A fundamental problem related to (1.1) is the relationship between p in O and the 
corresponding transmission eigenvalues. In [4,8] it has been observed that the transmission 
eigenvalues carry some information about p in O. The case n = 3 is naturally the most 
relevant in applications. A key question is whether we can uniquely determine p in O if 
all the transmission eigenvalues are known. Another important question is whether the 
unique recovery is possible if we know only a certain subset of the transmission eigenvalues. 

In the case where Q is the ball of radius b and p(x) is radially symmetric, it has recently 
been shown [9] that the set of all transmission eigenvalues uniquely determine p in O. In 
the radially symmetric case, let us use p{x) instead of p(x) with x := |x|. In this case it is 
natural to ask whether p(x) can be determined from a subset of transmission eigenvalues, 
such as those transmission eigenvalues for which the corresponding eigenfunctions are also 
spherically symmetric. We will refer to such eigenvalues as special transmission eigenvalues. 
Another variant of the transmission eigenvalue problem in the spherically-symmetric case 
has been studied in [24-26] , where some uniqueness results were established when only the 
positive special transmission eigenvalues are used in the determination. 

In the case n = 3, where O is the ball of radius b > and p is radially symmetric, the 
boundary- value problem (1.1) becomes equivalent to a nonstandard Sturm-Liouville-type 
eigenvalue problem, which is formulated in the following proposition. Here, "nonstandard" 
refers to the fact that the spectral parameter appears in the boundary condition at the right 
endpoint. Our assumptions on p are that p(x) is positive and continuously differentiable 
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and that p" is square integrable, i.e. 



p(x)>0, xe{0,b); peC\0,b); p" e £ 2 (0, b). 



(1.2) 



where a prime is used to denote the x-derivative. 

Proposition 1.1 Consider the special case of (1.1) with Vt being the three-dimensional 
ball of radius b centered at the origin, where only spherically- symmetric wave functions 
are allowed and it is assumed that such wave functions are continuous in the closure of 
O. Then, the corresponding special transmission eigenvalues of (1.1) coincide with the 
eigenvalues of the nonstandard boundary-value problem 

<&" + Xp(x) $ = 0, < x < 6, 

sin(v/A&) _ - 



d>(0) = 0, 



&(b)-cas(y/\b)$(b) = 0. 



PROOF: The Laplacian in R 3 in polar coordinates {x, 9, p) is given by 



1 9 2 d 
A := — - — x 2 — 



d . d 
smp— h 



d 2 



(1.4) 



x 2 dx dx x 2 sin cp dp dp x 2 sin 2 <p 09 2 ' 

where we recall that x := |x|. If the wave functions $ and $ are spherically symmetric, 
i.e. if they do not depend on 9 and p, then with the help of (1.4) we transform (1.1) into 



2*' 

\P" H + A p(x) * = 0, < x < 6, 



2*f 



(1.5) 



+ — - + A * = 0, < x < 6, 

I ¥(&) = ¥„(&), *'(&) = *&(&), 
where \I/ (0) and \I/o(0) must be finite because of the continuity of \& and \l/o i n ^- Letting 
$ := x^ and $ : = a^cb from (1.5) we get 

$" + Ap(x)$ = 0, 0<x<b, 

$% + \$ = 0, 0<x<6, (1.6) 
$(0) = $ (0) = 0, $(6) = $ (&), = $o(&)- 

4 



From the second line in (1.6) we see that the solution &o(x) satisfying <3>o(0) = must be 



The eigenvalues of (1.3), namely the A- values for which (1.3) has a nontrivial solution 
$(x), are the special transmission eigenvalues mentioned earlier. In other words, the 
corresponding eigenfunctions are spherically symmetric and hence functions of x only. Note 
that such eigenfunctions of (1.3) can only be determined up to a multiplicative constant, 
and it is clear from (1.3) that there exists only one linearly independent eigenfunction for 
each eigenvalue of (1.3). Nevertheless, for each eigenvalue \j of (1.3) we will associate a 
"multiplicity" in a special sense, namely the multiplicity of Xj as a zero of the quantity 
D(\) defined in (2.10). We will elaborate on the meaning of "multiplicity" in Section 2. 

We define the relevant quantity a as 



which has the physical interpretation as the travel time for the wave to move from x = 
to x = b. Our main result in this paper is the proof that the knowledge of eigenvalues of 
(1.3) with their "multiplicities" uniquely determine p(x) for < x < b provided a < b. If 
a = 6, we prove the unique determination of p provided we know one additional parameter, 
namely the value of the constant 7 appearing in (2.13). Let us clarify that we do not study 
the existence aspect of the inverse problem but we only analyze the uniqueness aspect. In 
other words, corresponding to our data we assume that there exists at least one function p 
satisfying (1.2), and we prove that if p\ and p^ are two such functions then we must have 
Pi = P2- 

When a = 6, it is an open question if knowledge of 7 is necessary or whether 7 can 
be determined from the knowledge of eigenvalues of (1.3) including their "multiplicities." 
In the discrete version of (1.3), assuming the existence aspect of the inverse problem is 
solved, it is already known [29] that generically, except for one exceptional case, p is 
uniquely determined from the knowledge of the special transmission eigenvalues and their 



a constant multiple of sin(v / A x)/\/~\. Thus, we see that (1.6) is equivalent to (1.3). | 




(1.7) 
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"multiplicities" and hence 7 is in general uniquely determined without needing to know 
any additional parameter. 

If p(x) satisfies (1.2), it is known [24,26] that 



where X nj for j e N are the real eigenvalues of (1.3) indexed in an increasing order, with 
N denoting the set of positive integers. Hence, the quantity a can be determined if the 
real eigenvalues of (1.3) are known. In other words, if p\ and P2 satisfy (1.2) and they 
correspond to the same set of special transmission eigenvalues, then we must have a\ = 0,2, 
where 



Let us elaborate on the eigenvalues of (1.3). As we illustrate with some examples in 
Section 2, besides real eigenvalues, (1.3) has in general nonreal eigenvalues and in fact the 
number of nonreal eigenvalues may be infinite. Because p(x) is real valued, from (1.3) it 
is seen that if A is an eigenvalue then A* is also an eigenvalue of (1.3), where we use an 
asterisk to denote complex conjugation. In our present work, for the unique recovery of p 
we assume the knowledge of all the eigenvalues (both real and complex nonreal) including 
their "multiplicities." In the previously established uniqueness results [24-26] regarding 
(1.3) it has been assumed that either a < 6/3 or that some partial information on p is 
available. On the other hand, in those results [24-26] it is assumed that only the positive 
eigenvalues are known and no "multiplicities" are used in the data. 

Our paper is organized as follows. In Section 2 we present some preliminary results 
that are needed to prove the uniqueness theorems of Sections 3 and 4. In Section 3 we 
consider the uniqueness in the recovery of p from the knowledge of special transmission 
eigenvalues of (1.1) with "multiplicities." When a < 6, where a is the quantity in (1.7), we 
establish the uniqueness. When a = 6, we show that the combined knowledge of special 



A 



(a - b) 2 



+ 0(1), 



Tij — > +OO, 




(1.8) 
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transmission eigenvalues of (1.1) with "multiplicities" and the constant 7 appearing in 
(2.13) assures the uniqueness. We also elaborate on the case a > b and indicate why the 
technique we use does not apply in that case to prove the uniqueness. In Section 4 we 
consider the uniqueness in the recovery of the potential V of the Schrodinger equation 
from the data consisting of special transmission eigenvalues of (4.1) with "multiplicities." 
We prove the unique recovery if our data contains one additional parameter, namely the 
constant 7 appearing in (4.5). 

2. PRELIMINARIES 

Let us recall [1] that an entire function of order 1/2 grows no faster than 0(e c l A l ( 1 ,+e ) 
as A — > 00 in the complex plane C for any given positive e, where c is some positive constant. 
The sums and products of such functions are entire of order not exceeding 1/2. 

We first consider a problem closely related to (1.3), namely 

<p" + Xp(x) = 0, 0<x<b, 

(2.1) 

0(0) = 0, 0'(O) = 1. 

It is known [30] that, for every A in the complex plane C, (2.1) has a unique solution 4>(x), 
which we also write as <f>(x; A) to emphasize its dependence on A. Since p(x) is real valued, 
the solution to (2.1) satisfies 

(j)(x;\*) = <f>(x;\)*, XeC. (2.2) 

Proposition 2.1 Assume that p satisfies (1.2). Then, (2.1) is uniquely solvable, and 
for each fixed x G (0,6] the quantities <f>(x; •) and <f>'(x; •) are entire in A of order 1/2. 
Furthermore, <f>(x; A) and <f>'(x; A) cannot simultaneously vanish at the same x-value. 

PROOF: We refer the reader to [30] for the proof that •) and <j>'(x; •) are entire in A of 
order 1/2. If 4>(xq] A) = 4>'(xo; A) = for some x$ value in [0, b], then the unique solution 
to the corresponding initial-value problem would have to be the zero solution, which is 
incompatible with <j)'(0) = 1 in (2.1). | 
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When p satisfies (1.2), it is known (see e.g. [30]) that the variable-speed wave equation 
in (2.1) can be transformed into a Schrodinger equation via a Liouville transformation. In 
other words, by using the change of variables 

y = y ( x ):= yfpjsjds, <p(y) = <p(y(x) ) : = p(x) 1/4 <j>{x) , (2.3) 
Jo 

we can transform (2.1) into the equivalent Sturm-Liouville problem for the Schrodinger 
equation that is given by 

-ip"(y) + q{y) ip(y) = A <p(y), < y < a, 



V (0) = 0, <p'(0) 



(2.4) 



where a is the quantity defined in (1.7) and 



1 p"(x) 5 p'{xf 



q(y) = q(y(x)) ■-- 



4p(x) 2 16 p{x) 3 



Let us use Im[\/A] to denote the imaginary part of \/A, where the argument of the 
square-root function is chosen so that arg(\/A) e (— n/2,n/2]. The proof of the following 
proposition can be obtained [30] with the help of the Liouville transformation (2.3) and 
some estimates for the corresponding Schrodinger equation in (2.4), and hence it will not 
be given here. 

Proposition 2.2 Assume that p satisfies (1.2). Then there exists a positive constant A 
such that, for all x G [0,6] and A G C, the solution <f>(x; A) to (2.1) and its x-derivative, 
respectively, satisfy 



0(x; A) 



[p(0) p(x)} 1/4 Vx 



sin 



(V\y(xf) <^exp(jlm[V\]\y(x)^ (2.5) 



\/A 



p{x) 
P(0) 



1/4 



COS 



V\y(x] 



< Aexp ^ Imfv/A] 



y(x] 



(2.6) 



where y(x) is the quantity given in (2.3). 



For a positive e let C e denote the sector in the complex plane defined as 

C e := {A G C : e < arg(A) < 2tt - s}. 



(2.7) 



The proof of the following result is already known [27]. 

Proposition 2.3 Assume that p satisfies (1.2). Then, for each fixed x G [0, 6], as A — > oo 

in C g , the unique solution 4>(x; A) to ^.ij satisfies 



<f>(x; A) 



[p(0)pOr)] 1/4 VA 



sin 



v / Ay( 



1 + 0| 7x 



/(*;A) = 



>(x)" 


1/4 


Lp(o). 





COS 



y/Xy{x)^j 



1 + 



v/A 



(2.8) 
(2.9) 



where y(x) is the quantity given in (2.3). 

Let us now clarify the relationship between (1.3) and (2.1). In general, for a given 
A G C, (1.3) may not have a nontrivial solution. Suppose that Xj is an eigenvalue of (1.3). 
Then a solution 3>(x; Xj) to (1.3) can only be determined up to a multiplicative constant, 
and in fact any such solution must be a constant multiple of the unique solution <f>(x; Xj) 
to (2.1) due to the fact that $(0) = in (1.3) and 0(0) = in (2.1). 

We now introduce the key function D(X) as 

sin(v^A6) . 



D(X) 



y/X 



<f>'(b; A) - cos(a/A6) 4>{b] A), 



(2.10) 



where we recall that <j>(x; A) is the unique solution to (2.1). Let us remark that, if p{x) = 1 
in (2.1), then <f)(x; A) = sin(v^A x)/VX and hence D(X) = 0. 

Theorem 2.4 Assume that p satisfies (1.2). Then, the quantity D(X) defined in (2.10) 
is entire in X of order not exceeding 1/2. Each zero of D(X) in the complex plane C 
corresponds to an eigenvalue of (1.3) and vice versa. The value X = is always a zero of 
D(X) of some multiplicity d with d > 1, and hence 



D(0) = 0. 
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(2.11) 



Furthermore, 



D(X*) =D(X)*, 



A G C, 



(2.12) 



and there exists a real constant 7 such that 



D(A) = 7 S(A), 



(2.13) 



where the auxiliary quantity S(A) is uniquely determined from the zeros (including multi- 
plicities) of D(X) and has the representation 



with X n for n G N being the nonzero zeros of D(X), some of which may be repeated. 
PROOF: From their representations in terms of exponential functions, we know that 



that 4>(b; A) and (j>'(b; A) are entire of order 1/2, and hence the right side of (2.10) is entire 
of order not exceeding 1/2. If Xj is an eigenvalue of (1.3) with an eigenfunction 3>(x; Xj), 
we already know that $(x; Xj) is a constant multiple of the solution (f>(x; Xj) to (2.1), and 
hence from (2.10) we see that D(Xj) = 0. Conversely, if D(Xj) = for some Xj, then a 
comparison of (1.3) and (2.1) shows that the unique solution (f)(x; Xj) to (2.1) satisfies (1.3) 
and hence Xj is an eigenvalue for (1.3) with eigenfunction <f>(x; Xj). In particular, we note 
that when A = the unique solution to (2.1) is given by 




(2.14) 



sin (VXb)/VX and cos(v / A6) are entire in A of order 1/2. From Proposition 2.1 we know 



(f)(x; 0) = x, 



(2.15) 



which indicates that 



0(O;O) = O, 0'(O;O) = 1, <£(&; 0) = 6, 0'(&;O) = 1, 



(2.16) 



and hence (f>(x; 0) indeed satisfies (1.3) when A = 0. Thus, A = is always a zero of D(X) 
with some multiplicity d, which is at least one. We obtain (2.12) from (2.2) and (2.10). 

10 



Since D(X) is entire of order not exceeding 1/2, by the Hadamard factorization theorem, 
we must have the representation in (2.13), where 7 is a complex constant and S(A) as in 
(2.14). In fact 7 turns out to be real as a result of (2.12). | 

As we have seen in Proposition 1.1 and Theorem 2.4, the special transmission eigen- 
values of (1.1), the eigenvalues of (1.3), and the zeros of D(X) defined in (2.10) all coincide. 
On the other hand, each zero of D(X) may have a multiplicity greater than one even though 
there exists only one linearly independent eigenfunction for the corresponding eigenvalue 
of (1.3). We refer to the multiplicity of a zero Xj of D(X) also as the "multiplicity" of 
the special transmission eigenvalue Xj. Next, we elaborate on the "multiplicities" with an 
illustrative example. 

Example 2.5 When p(x) is constant on (0,6), by using p to denote that constant value, 
the unique solution to (1.3) is obtained as 

<f)(x; X) = —= sm(\/~\px), < x < 6, 
V-V 

and hence the corresponding quantity in (2.10) is given by 

1 1 

D(X) = —= sin(v / A6) cos( v / V&) - -7= cos(v / A6) sm(^/Xpb). (2.17) 
VA y/Xp 



When p(x) = 1/4, from (2.17) we get 



2 . 3 ( VXb 



D(X) = ^ sin 

and hence D(X) has a simple zero at A = and an infinite set of real zeros at the A- 
values Aj 2 n 2 /b 2 for j e N, each having a multiplicity of three. On the other hand, when 
p(x) = 4/9, from (2.17) we get 



D(X) = —= sin 



1 .of VXb 



3 + 2 cos 



2VXb 



VX \ 3 

and hence D(X) has a simple zero at A = 0, an infinite set of real zeros of multiplicity 
three at the A-values 9j 2 7r 2 /6 2 for j e N, and an infinite set of simple complex zeros at 
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the A- values that are given by 



9(2j - l) 2 n 2 9 



Ab 2 



Ab 2 



log 



±i 



9(2j - 1)tt 
2b 2 



log 



'3 + \/5 X 



j GN. 



Let us remark that the knowledge of S(A) given in (2.15) is equivalent to the knowledge 
of the eigenvalues of (1.3) with their "multiplicities." Furthermore, the knowledge of H(A) is 
equivalent to the knowledge of its zeros including their multiplicities. Hence, in proving our 
uniqueness results, as our data we can equivalently use H(A), the zeros of H(A) with their 
multiplicities, the eigenvalues of (1.3) with their "multiplicities," or the special transmission 
eigenvalues of (1.1) with their "multiplicities." 

Since D(X) given in (2.10) is entire, we can obtain its Maclaurin expansion with the 
help of the Maclaurin expansion of the unique solution A) to (2.1), which we write as 



4>(x\ A) = <f> (x) + A0i(x) + \ 2 (f)2{x) + 0(A 2 ), A ->■ in C, 



(2.18) 



where we have defined 



4> (x) := x, 4>i(x) := M 2 (x) - xM^x), 

2 rX t>x r>x 

(f) 2 (x):=- dziM^z)] 2 -x / dzzM 2 (z)+ dz z p(z) M 2 (z), 
2 Jo Jo Jo 



with 



Mi(x) := / dzzp(z), M 2 (x) := / dzz 2 p(z). 
Jo Jo 

Using (2.18)-(2.20) and their ^-derivatives in (2.10) we obtain 



(2.19) 
(2.20) 



D{\) =D + XD 1 + X 2 D 2 + 0(X S 



A ->■ in C, 



(2.21) 



where 



Do 



D o :=0, D 1 :=j-M 2 (b), 
~+b[M 1 (b)] 2 -M l (b)M 2 (b)-jM 1 (b) + b ^M 2 (b)- I dzlM^z)]' 
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If d = 1 in (2.14), with the help of (2.10) and (2.21) we see that 

L»!^0, 7 = £>i, - 7 V— = J D 2 , 

where 7 is the parameter appearing in (2.13) and Xj for j G N are the nonzero zeros of 
-D(A), some of which may be repeated. On the other hand, if d = 2 in (2.14), then we must 
have 

D 1 = 0, D 2 ^ 0, 7 = £> 2 . 

The results in the following propositions will be used in the proof of the unique 
determination of p. 

Proposition 2.6 Suppose that h is an entire function of A, and let 

/(A):=! ta^) h(A)j (222) 

for some positive constant c, and assume that /(A) is entire of order not exceeding 1/2. 
Then, the order of h(X) cannot exceed 1/2. 

PROOF: Note that (2.22) and the fact that the order of / does not exceed 1/2 imply that 
for any positive e there exists a positive constant A such that 

H\)\ < , r. I/(A)| < , . T~7t r, U\X\^) . (2.23) 

|sm(vAc)| |sm(vAc)| v J 

In the neighborhood of the zeros of sin(v / Ac)/v / A, which occur when A = n 2 n 2 /c 2 for 
n G N, the bound in (2.23) is too large to assure that the order of h(X) cannot exceed 1/2. 
Thus, we need to analyze the behavior of h(X) near those zeros. Let us enclose each such 
zero within the disk U n of radius one, where we have defined 



U n :=<\eC 



2 2 

X n— 

c z 



< 1\ , n e N. 



The boundary dll n can be parameterized by using 1} so that if A G dU n then 



X=^ + e i », ^G(-7r,7r], 
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or equivalently 



From (2.24) we get 



c 2 e^ 1 



\c z = n TV ( 1 + — =- , AG 8U n . 



K c 2 e^ 1 / 1 

VXc = nn+ — + O — , 

In n \n ,i 



A G 9C/ n , n +oo, 



where we recall that arg(A) G (— 7r/2,7r/2]. Using the trigonometric relations 



sin(z + rar) = (-l) n sin(z), 2 G C, 



(2.24) 



(2.25) 



sin(z) = z + o(z), z — > in C, 



from (2.25) we get 



sin(\/A ( 



c z 1 /l 
— - + o - 

2n n \n 



A G <9L/n, n — >■ +00. 



Note that (2.25) implies that 

VXc 



n = 



A G dU n , A — >■ 00, 



and hence we can write (2.26) as 



sin(v^A c) 



+ o(-j=), A G dU n , A ->■ 00. 



(2.26) 



(2.27) 



21^1 Vv/A, 

The estimate in (2.27), the fact that sin(\/Ac) grows exponentially for large Im[\/A] while 
its nonzero zeros are confined to the centers of the disks U n , and the minimum modulus 
principle applied to the exterior of U n imply that there exist positive constants m and M 
such that 



sin(v / A < 



> 



m 



|Va| 



agc\u»a, |a|>m, 



(2.28) 



and hence (2.23) and (2.28) yield 
A\V\\ 



\h{\)\ < 



m 



exp (c|A|( 1/2)+e ) , \eC\U™ =1 U n , \\\>M. 



(2.29) 
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On the other hand, by the maximum modulus principle, the maximum of \h(X)\ in the 
closure of U n must occur on the boundary dU n , and hence (2.29) holds, whenever |A| > M, 
perhaps by replacing A there with another positive constant. Hence, we have proved that 
the order of h cannot exceed 1/2. | 

Proposition 2.7 Let f be an entire function of X such that 



< 




A — > oo in C, 



(2.30) 



where c is a positive constant. Then there is a constant C\ such that 




(2.31) 



Similarly, if g is an entire function of A such that 




(2.32) 



then there is there is a constant C2 such that 




(2.33) 



PROOF: The second line of (2.30) implies that /(A) can be written as 




(2.34) 



for some entire function h(X). Using (2.34) in the first line of (2.30) we get 




A G C, 



(2.35) 
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for some positive constant B. By Proposition 2.6 we know that the order of h cannot exceed 
1/2. On the other hand, by using the exponential representation of the sine function, as 
A — ?■ oo along any ray other than the positive real axis we have 

sin ( y/\c 



Vx 



exp ^ 


Im[VX] 




2 


v/A 





[1 + 0(1)]. 



(2.36) 



Hence, from (2.35) and (2.36) we see that h(X) must be bounded on any ray other than 
the positive real axis. By invoking a consequence of the Phragmen-Lindelof principle (see 
Theorem 18.1.3 of [17]) we conclude that h(X) must be a constant, which establishes (2.31). 
The proof of (2.33) is obtained in a similar manner. | 

Proposition 2.8 Let f be an entire function of X satisfying (2.30), and assume that as 
X — > oo along some fixed ray in the complex plane we have 

exp ( Im[\/A] 



/(A) = 



Vx 



o(l). 



(2.37) 



Then, /(A) = 0. Similarly, let g be an entire function of X satisfying (2.32), and assume 
that as X — > oo along some fixed ray in the complex plane we have 



g(X) = exp ( Im[VX] cj o(l) 



(2.38) 



Then, g(X) = 0. 

PROOF: In the proof of Proposition 2.7, the further restriction given in (2.37) forces us to 
have C\ = in (2.31), and hence we get /(A) = 0. Similarly, (2.38) forces to have C2 = 
in (2.33), yielding g{X) = 0. | 

We state a relevant relationship between (2.1) and two Sturm-Liouville problems in 
the following corollary. 

Corollary 2.9 Let p satisfy (1.2). Then, the eigenvalues of the Sturm-Liouville problem 

ip" + X p(x) i) = 0, < x < b, 

^(o) = m = 0, 

16 



(2.39) 



exactly correspond to the zeros of (/>(&; A), where <f>(x; A) is i/ie unique solution to (2.1). 
Similarly, the eigenvalues of the Sturm- Liouville problem 



exactly correspond to the zeros of 4>'(b; A). 

The fundamental uniqueness theorem of inverse spectral theory for Sturm-Liouville 
problems indicates that, assuming the existence problem is solved, the knowledge of two 
sets of spectra uniquely determines p. It is already known [2,16,19-23] that the combined 
knowledge of the eigenvalues of (2.39) and the eigenvalues of (2.40) uniquely determines 
p{x) for x G [0, b]. Thus, with the help of Corollary 2.9 we have the following result. 

Corollary 2.10 Let pi and p^ satisfy (1.2), and let 4>i(x; X) and foix; A), respectively, be 
the corresponding unique solutions to (2.1). Then, p\ = p2 if <j>i(b; X) and ^2(6; A) have 
the same set of zeros and also (^[(b; A) and (f) 2 (b', A) have the same set of zeros. 

3. THE INVERSE PROBLEM 

We assume that p satisfies (1.2). The relevant direct problem is the determination 
of the special transmission eigenvalues of (1.1) including their "multiplicities" when p(x) 
is known for x G [0,6]. Conversely, our relevant inverse problem is the determination of 
p(x) for x G [0,6] from the knowledge of the special transmission eigenvalues of (1.1) 
including their "multiplicities." From (2.13) and (2.14) we see that the direct problem can 
be equivalently stated as the determination of the map p (->■ H and the inverse problem as 
the determination of the map S h-> p, where S is the quantity appearing in (2.14). Recall 
that we are only concerned with the uniqueness aspect of the inverse problem and not 
with the existence aspect. In other words, corresponding to our data we assume that there 
exists at least one function p satisfying (1.2) and we show that our data leads to a unique 
p. 




< x < 6, 



(2.40) 
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The main conclusion in our paper is that, once the existence problem is known to be 
solvable, the function S uniquely determines p in case a < 6, where a is the quantity defined 
in (1.7). On the other hand, when a = 6 it is unclear if S uniquely determines p, but we 
show that S and 7 together uniquely determine p, where 7 is the constant appearing in 
(2.13). In other words, if a = b then p(x) for x G [0, b] is uniquely determined by -D(A) for 
A G C. First, we present a special case of the uniqueness result in the following theorem, 
which also includes the solution to the relevant existence problem. 

Theorem 3.1 Assume that p satisfies (1.2), andletthe corresponding D{\) be as in (2.10). 
IfD(X) = for A G C, then p(x) = 1 for x G [0,6]. 

PROOF: If D(X) = then (2.10) implies that 

^^2^0'(6;A) = cos(v / A6)0(6;A), A G C. (3.1) 
v A 

Note that each of the four functions in (3.1), namely, sin(-\/A b)/y/\, cos(v^A6), 0(6; A), and 
0'(6; A) are entire of order 1/2. Furthermore, 0(6; A) and 0'(6; A) cannot vanish simultane- 
ously, and sin(v^A6)/-\/A and cos(\/A6) cannot vanish simultaneously. Thus, (3.1) implies 
that sin(v^A6)/v^A and 0(6; A) must have the same set of zeros including multiplicities and 
that cos(-\/A6) and 0'(6; A) must have the same set of zeros including multiplicities (note 
that, in this particular case, the multiplicities must all be one). Hence, by the Hadamard 
factorization theorem, considering the fact that the order of each of these four functions is 
1/2, we must have 

0(6; A) = ci sm (^ b ) ; 0'(6;A) =ci cos(v / A6), (3.2) 
v A 

for some nonzero constant C\\ in fact, (2.16) implies that C\ = 1. By Corollary 2.10 we know 
that the combined knowledge of the zeros of 0(6; A) and of 0'(6; A) uniquely determines p. 
Thus, p is uniquely determined by the combined knowledge of the zeros of sin(-\/A b)/\f\ 
and of cos(\/A 6), and it is already known that those combined zeros correspond to p{x) = 1 
for x G [0,6]. I 
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In the next theorem, we present our uniqueness result when a < b. 



Theorem 3.2 Assume that for the function S appearing in (2.14) there corresponds at 
least one function p satisfying (1.2); assume also that a < 6, where a is the quantity defined 
in (1.7). Then, p is uniquely determined by S; in other words, the knowledge of special 
transmission eigenvalues of (1.1) with "multiplicities" uniquely determines p. 

PROOF: Let us assume that p\ and p 2 correspond to Si and S 2 , respectively, and let D 1 (X) 
and D 2 (\) be the corresponding quantities in (2.13) with 71 and 72 being the respective 
constants there. We will show that p\ = p 2 if Si = S 2 . Let us also use (pi and (p 2 to 
denote the solutions to (2.1) corresponding to p\ and p 2 , respectively. From the line above 
(1.8) we see that if Si (A) = S 2 (A) then a\ = a 2 , where a± and a 2 are the corresponding 
quantities for p\ and p 2 , respectively. Let us use o to denote the common value of a\ and 
a 2 . Since we assume that a < 6, by (2.5), (2.6), (2.8), and (2.9) we have 

exp 



ii(6; A) = 



>i(6;A) = 



exp 



v/A 
(| Mv^A] 



a ) exp 

0(1), 2 (6;A) = V 

) exp ( \hn[V\] 

o(l), ( p 2 (b;X) = — 



Vx 



0(1), A^ooinC, 
(3.3) 

o(l), A — ?■ 00 in C e , 
(3.4) 

0'i(6; A) = exp ( Im[y/\] a) 0(1), <p 2 (b; A) = exp ( Im[y/\] a) 0(1), A ->■ 00 in C, 

(3.5) 

</>;(&; A) = exp ( M-v/X] 6 ) <^( & ; A ) = ex P ( Imfv/X] ft) o(l), A ->■ 00 in C e , 

(3.6) 

where C e is the sector defined in (2.7). Since Si (A) = S 2 (A), from (2.13) and (2.14) it 
follows that 



2 2 



and hence from (2.10) we get 

1 



2 2 



1 A Ih nV 



neN. 



(3.8) 
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In a similar way, with the help of (2.9), (2.10), (2.13), and (2.14), by using 




) 



neN, 



(3.9) 



we obtain 




(2n-l) 2 7T 2 
' 4b 2 



) 



— 02 ( b ; 

72 V 



(2n-l) 2 7T 2 
' 4b 2 



) 



n G N. 



(3.10) 



In Proposition 2.8 by choosing c = b and /(A) = 4>i(b; \)/"fi — 4>2(b; A)/72, we see (3.3), 
(3.4), and (3.8) imply that /(A) = 0. Similarly, in Proposition 2.8 by choosing c = b and 
g(X) = (f)'i(b; A)/7i — 4>' 2 (b; A)/72, we see that (3.5), (3.6), and (3.10) imply that #(A) = 0. 
On the other hand, /(A) = indicates that 4>i{b; A) and <j>2(b; A) have the same set of zeros, 
and g(X) = indicates that 01(6; A) and </> 2 (6; A) have the same set of zeros. Thus, using 
Corollary 2.10 we conclude that p\ = p 2 - I 

The next uniqueness theorem applies to the case a = b. 

Theorem 3.3 Assume that for the function S appearing in (2.14) there corresponds at 
least one function p satisfying (1.2); assume also that a = 6, where a is the quantity defined 
in (1.7). Then, p is uniquely determined by the combined knowledge o/H and the constant 
7 appearing in (2.13); in other words, the knowledge of special transmission eigenvalues of 
(1.1) with "multiplicities" along with the knowledge 0/7 uniquely determines p. 

PROOF: The proof is similar to the proof of Theorem 3.2 with appropriate modifications 
we indicate here. As in the proof of Theorem 3.2 we have Si (A) = S 2 (A), but we also have 
7i = 72, and we want to show that p\ = pi- By Corollary 2.10 it is sufficient to prove that 
4>i (6; A) and 2 (6;A) have the same set of zeros and that 0^(6; A) and 2 (6;A) have the 
same set of zeros. Since a = 6, this time we have (3.3) and (3.5), but not (3.4) or (3.6). 
Proceeding as in (3.7)-(3.10) verbatim, and in Proposition 2.8 by choosing / and g as in 
the proof of Theorem 3.2, we obtain 



-Mb; A) -—Mb; A) = Ci 




(3.11) 




x/A 
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— (b; A) - — 0a(6; A) = C 2 cos(>/A&), 
7i 72 



(3.12) 



for some constants C\ and C 2 . Evaluating (3.11) and (3.12) at A = and using (2.16), we 
get 

d = C 2 = — - — . (3.13) 

7i 72 

Since we assume 71 = 72, we see from (3.13) that C\ = C2 = 0. Thus, from (3.11) 
and (3.12) we get 0i(6; A) = 02 (fy A) and 0^(6; A) = (f>' 2 (b; A), indicating that 0i(6; A) and 
02(6; A) have the same set of zeros and also that 0i(6; A) and (j)' 2 {b; A) have the same set of 
zeros. I 

Having considered the inverse problem when a < b and a = b in Theorems 3.2 and 
3.3, respectively, let us now comment on the case a > b. The method we use to prove the 
uniqueness for a < b does not apply to the case a > 6, as the following argument clarifies. 
The lack of applicability of our technique to the case a > b certainly does not mean that a 
uniqueness result does not exist when a > b. The unique recovery of p from D(X) defined 
in (2.10) is based on our ability to extract each of 4>(b; X) and </>'(&; A) up to a constant 
multiplicative factor. When a > b let us the consider the determination of two functions 
4>(b; X) and <j>'(b;X) that are entire in A and of order 1/2 and that satisfy the respective 
asymptotics related to (2.5) and (2.6), namely, as A — > 00 in C 



<P(b; A) 
for which 



exp ( Im[\/A] a J , , 
L 0(1), cf>'(b; A) = exp ( Im[VX] a) 0(1), 



Vx 

sin(v / A6) 



(b; A) - cos(v / A6) <p(b; A) = D(X). 



Let C be any entire function of A having the asymptotics 



C(A) = exp ( Im[\^X] (a - 6)) 0(1), A ->• 00 in C. 



(3.14) 



(3.15) 



(3.16) 



Letting 



<P(b; A) :=0(6; A) + 



S111< V A h) C(A), A) := <p'(b; A) + cos(v / A6) C(A), 
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we see that </>(&; A) and <//(&; A) are entire in A and that (3.14) and (3.15) are satisfied 
when we replace in them 0(6; A) with 0(6; A) and replace <j>'(b; A) with </>'(&; A). Because of 
(3.16), £(A) must be a constant when a = b and must be zero when a < 6, but no such 
restrictions exist when a > b. Thus, our method does not allow us to conclude the unique 
determination of p from D(X) when a > b. 

4. THE INVERSE PROBLEM FOR THE SCHRODINGER EQUATION 

In the case of the Schrodinger operator, the interior transmission eigenvalue problem is 
analogous to the corresponding problem for the wave equation with variable speed. Instead 
of (1.1), we have 

-A* + V(x) * = x G O, 

on on 

where /U is the spectral parameter, V(x) is a real- valued potential that is square integrable 
on O, and it is assumed that V(x) = outside O. Those /i-values yielding nontrivial 
solutions \I/ and \I/o to (4.1) are called transmission eigenvalues of (4.1). In the spherically- 
symmetric case, using V(x) instead of V(x) with x := |x|, we have the following analog of 
Proposition 1.1. We omit its proof because it is similar to the proof of Proposition 1.1. 

Proposition 4.1 Consider the special case of (4-1) with being the three-dimensional 
ball of radius b centered at the origin, where only spherically- symmetric wave functions 
are allowed and it is assumed that such wave functions are continuous in the closure of O. 
Then, the corresponding transmission eigenvalues of (4-1) coincide with the eigenvalues of 
the boundary-value problem 

+ V(x) $ = <x < 6, 

$( ) = 0, Sln(v/M) $'(&) - cos( v ^6) $(&) = 0. (4 ' 2) 

The eigenvalues of (4.2), namely the //-values for which (4.2) has a nontrivial solu- 
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tion coincide with the special transmission eigenvalues of (4.1), namely those transmission 
eigenvalues of (4.1) for which the corresponding wavefunctions are spherically symmetric 
in addition to V being spherically symmetric. Note that the boundary condition at x = b 
in (4.2) suggests an analog of D(\) appearing in (2.10). We define 

sin( v //I6) ~ 



VP 



<p'(b; fj) - cos(y/jib)<i>(b;n), 



(4.3) 



where <f>(x; Li) is the analog of <f>(x; A) appearing in (2.4) and is the unique solution to the 

initial-value problem _ _ 

" -(/)" + V(x) (j) = ii (j), <x <b, 

(4-4) 

0(0) = 0, 0'(O) = 1. 
The following proposition contains results that are analogous to those stated in Propo- 
sitions 2.1, 2.2, and 2.3. 

Proposition 4.2 Assume that V is real valued and square integrable on [0,6]. Then, for 
each /i G C, (4-4) has a unique solution (f)(x; y). For each fixed fj,, the functions </>(•; y) and 
</>'(•; (j,) cannot vanish at the same x-value. For each fixed x G (0, b], the quantities <j>(x; •) 
and (f)'(x; •) are entire in fi and there exists a positive constant A such that 



4>(x;ii) ^— sin(y/JIx) 

v t 1 



< 



A 



Iv^l 



expdlmf^l x) , 



'{x; n) — cos {y/Jix) < Aexp (\Im[y/Jj]\ x) . 

Furthermore, for each fixed x e [0, 6], as /i — > oo in the sector C e defined in (2.7), we have 

1 



<f>(x;n) = — sin (y/jix) 



1 + 



4>'(x; fl) = cos (y/Jix) 



l + O 



1 

v 7 ^ 



We note that (4.2), (4.3), and (4.4) are closely related. If fj,j is an eigenvalue of (4.2) 
with an eigenfunction 3>(x;iij), then 3>(x; lij) must be a constant multiple of (f)(x;/ij), 
where fi) denotes the unique solution to (4.4). Hence, from (4.3) we conclude that 
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D(fj,j) = 0. Thus, with the help of Proposition 4.1 we conclude that the special transmission 
eigenvalues of (4.1), the eigenvalues of (4.2), and the zeros of (4.4) all coincide. On the 
other hand, there exists only one linearly independent eigenfunction for a given eigenvalue 
Hj of (4.2) whereas the multiplicity of fij as a zero of D(/j) may be greater than one. 
We will refer to the multiplicity of a zero fij of D{fi) as the "multiplicity" of the special 
transmission eigenvalue fij and also as the "multiplicity" of the eigenvalue fij of (4.2). 

Note that from (4.3) and the second line of (4.4) we obtain 

£(0) = &#'(&; 0)-#(6;0), 

and, contrary to (2.11), generically we have D(0) ^ 0, although we may have -D(O) = for 
some potentials. For example, if V(x) = 0, then we have (f)(x; fi) = sm(^JIx)/ \JJi, yielding 
-D(O) = 0. In fact, D(fi) = for fi e C in that special case. 

Our goal in this section is to show that V(x) for < x < b is uniquely determined by 
the corresponding D(fi) known for all fi e C. In fact, we will see that, up to the multiplica- 
tive constant 7 appearing in (4.5), D(fi) is uniquely determined by the knowledge of its 
zeros including the multiplicities of those zeros. Since those zeros are exactly the eigenval- 
ues of (4.2), we will conclude that the knowledge of the eigenvalues of (4.2) including their 
"multiplicities" and the value of 7 uniquely determines V. Since the eigenvalues of (4.2) 
are the special transmission eigenvalues of (4.1), we will also conclude that the knowledge 
of those special transmission eigenvalues including their "multiplicities" and 7 uniquely 
determines V. Since the proofs are similar to those in the case of the variable-speed wave 
equation studied in the previous sections, we will omit some of the proofs. 

As in the case of the variable-speed wave equation, we only consider the uniqueness 
aspect of our inverse problem and not the existence aspect. In other words, corresponding 
to our data D(fi) or its equivalents, we assume that there exists at least one potential V, 
where V(x) is real valued and belongs to L 2 (0, b). We then prove that if Vi and V2 are two 
such potentials, then we must have V\ = V2. Let us also clarify that the equality V\ = V2 is 
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meant to be an equality in the almost-everywhere sense because we deal with potentials in 
the class L 2 (0,6) whereas the corresponding equality p\ = P2 obtained in Section 3 holds 
pointwise because p\ and pi satisfy (1.2). 

The following theorem summarizes the properties of D{p) defined in (4.3), and it is an 
analog of Theorem 2.4. We omit the proof because it is similar to the proof of Theorem 2.4. 

Proposition 4.3 Assume that V is real valued and square integrable on [0,6]. Then, the 
quantity D(p) defined in (4-3) is entire in p and its order does not exceed 1/2. Thus, by 
the Hadamard factorization theorem D(p) is determined, uniquely up to a multiplicative 
constant, from its zeros as 



oo 



D(p)=7P d H 1- — ), (4-5) 

n=l V ^ nJ 



with p n for rgN being the nonzero zeros of D(p,), some of which may be repeated, and d 
denoting the multiplicity of the zero as a zero of D(p). 

The results stated in the following theorem are analogous to those stated in Corollar- 
ies 2.9 and 2.10. 

Theorem 4.4 Assume that V is real valued and square integrable on [0, 6], and let <ft(x, p) 
denote the unique solution to (4-4)- We then have the following: 

(i) The zeros of 4>(b; p) coincide with the eigenvalues of the Sturm- Liouville problem 

-i)" + V{x) V> = p, < x < b, 

(4-6) 

^(0) =tp{b) =0. 

(ii) The zeros of(f>'(b;p) coincide with the eigenvalues of the Sturm- Liouville problem 

+ V(x) V> = p$, 0<x<b, 

(4-7) 

^,(0) = V'(6) =0. 

(iii) The data consisting of the eigenvalues of (4-6) and (4-7) uniquely determines V if 
the existence is assured. In other words, assuming that there exists at least one V 
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corresponding to the data, if V\ and V2 correspond to the same data then we must 
have Vi(x) = V2(x) on [0,6]. 

(iv) The data consisting of the zeros of 0(6; y) and <fi'(b;fi) uniquely determine V if the 
existence is assured. In other words, assuming that there exists at least one V cor- 
responding to the data, if V\ and V2 correspond to the same data then we must have 
Vi(a;) = V 2 (x) on [0,6]. 

PROOF: We obtain (i) and (ii) by comparing (4.4) and (4.6) and by noting that 0(0) = 
ip(0) = 0. We note that (iii) is a version of the well-known uniqueness result by Borg [3]. 
Finally, (iv) is a consequence of (i)-(iii). | 

The following is the analog of the uniqueness result stated in Theorem 3.3. 

Theorem 4.5 Assume that V(x) is real valued and square integrable on [0,6]. Then, V is 
uniquely determined by the function D(/i) appearing in (4-3) if we assume that there exists 
at least one V corresponding to D. Equivalently stated, if the existence is assured, V is 
uniquely determined by the knowledge of the special transmission eigenvalues of (4-1) with 
their "multiplicities" and the constant^/ appearing in (4-5). 

PROOF: If Vi and V2 correspond to -Di(^) and D 2 (n), then we need to show that V\ = V2 
when Di{ji) = .D2 (/-*)• Let (fii(x; fi) and (f>2(x; ji) be the solutions to (4.4) corresponding 
to Vi and V2, respectively. Because of Theorem 4.4 (iv), it is sufficient to show that 
01 (6; (j,) = 02(6; y) and 0^(6; y) = 2 (6; fi), which is proved by proceeding as in the proof 
of Theorem 3.3. | 

By Proposition 4.3 we know that the knowledge of D(n) is equivalent to the knowledge 
of its zeros with their multiplicities and the constant 7 appearing in (4.5). We have 
already seen that the zeros of D{p), the eigenvalues of (4.2), and the special transmission 
eigenvalues of (4.1) all coincide. Thus, from Theorem 4.5 we obtain the following corollary. 

Corollary 4.6 Assume that V(x) is real valued and square integrable on [0,6]. Assuming 
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that there exists at least one V corresponding to the data, V is uniquely determined by the 
data consisting of the zeros of D(fi) in (4-5) with their multiplicities and the constant 7 
there. Equivalently, assuming the existence, V is uniquely determined by the data consisting 
of the eigenvalues of (4-%) with their "multiplicities" and the constant 7. 

One consequence of Corollary 4.6 is that if D(pt) = 0, then V(x) = 0, which is the 
analog of Theorem 3.1. 

Let us mention that it is an open problem whether the value of 7 appearing in (4.5) 
can be determined from the zeros of D{ji). If the answer is yes, then 7 is not needed for 
the unique determination of V, and as seen from Corollary 4.6 the zeros of D(fj) with 
their multiplicities would be sufficient for that purpose. The technique we use to prove the 
uniqueness assumes the knowledge of 7, but this does not rule out the possibility that there 
might be another method to obtain the uniqueness from the data consisting only of the 
zeros of D(fi) and their multiplicities. We note that in the discrete version of the inverse 
transmission problem for the Schrddinger equation, 7 is determined [29] in the generic case 
by the zeros of D(y) and their multiplicities. 
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